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Abstract According to the character of Gaussian, we 
modify the Rank-Ordered Absolute Differences (ROAD) 
to Rank-Ordered Absolute Differences of mixture of 
Gaussian and impulse noises (ROADG). It will be more 
effective to detect impulse noise when the impulse is 
mixed with Gaussian noise. Combining rightly the ROADG 
with Optimal Weights Filter (OWF), we obtain a new 
method to deal with the mixed noise, called Optimal 
Weights Mixed Filter (OWMF). The simulation results 
show that the method is effective to remove the mixed 
noise. 
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1 Introduction 

Noise can be systematically introduced into digitized 
images during acquisition and transmission, which usu- 
ally degrade the quality of digitized images. However, 
various image-related applications, such as aerospace, 
medical image analysis, object detection etc., generally 
require effective noise suppression to produce reliable 
results. The problem of noise removal from a digitized 
image is one of the most important ones in digital image 
processing. The nature of the problem depends on the 
type of noise to the image. Generally, two noise models 
can adequately represent most noise added to images. 
Often in practice it is assumed that the noise has two 
components: an additive Gaussian noise and an impulse 
noise. 

The additive Gaussian noise model is: 



Y(x) = f(x) + e(x), xel, 



(1) 



where I 



,^,1} 2 , N G N, Y is the 
observed image brightness, / : I h+ [c, d\ is an un- 
known target regression function, and e(x), x G I, are 
independent and identically distributed (i.i.d.) Gaus- 
sian random variables with mean and standard de- 
viation a > 0. The additive Gaussian noise is charac- 
terized by adding to each digitized image pixel a value 
from a zero-mean Gaussian distribution. Such noise is 
usually introduced during image acquisition. The zero- 
mean property of this Gaussian distribution makes it 
possible to remove the Gaussian noise by Non-Local 
weighted averaging. Important denoising methods for 
the Gaussian noise model have been well developed 
in rec e nt years, see for example Pol zchl an d Spokoinv 
(|200r | iBuades et al.l (120051) 



Kervrann and Boulanger 
20061) . I Aharon et al.l(|2006l ). lDabov et all (|2007h . lCai et al 
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(2008). Hammond and Simoncelli (2008 


). Lou et al. 


Katkovnik et al.l (hoich. and 1 Jin et al. 
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The random impulse noise model is: for x G I, 



Y{x) 



n{x), if x G B, 
/(x),ifxGl\B, 



(2) 



1996), 



(2004), 



for all pixels in the image, which take value in the inter- 
val (0, 1]. The weight will get low value (the value may 
be near to 0) when a pixel is contaminated by impulse 
noise; otherwise, the weight will carry a high value (the 
value may be near to 1). Then the Optimal Weights Fil- 
ter (OWF) combining with the ROADG statistic can 
detect the impulse points in the image and give the 
proper weights to deal with the mixed noise. As a re- 
sult, we obtain our new filter. The simulation results 
show that the proposed filter can effectively remove the 
mixture of impulse noise and the Gaussian noise. More- 
over, when applied to either the single impulse noise or 
the single Gaussian noise it performs as good as the 
best filters specialized to single noises. 

The rest of the paper is organized as follows. In Sec- 
tion[2]after a short recall of the Optimal Weights Filter 
and a brief presentation of the Trilateral Filter whose 
.Abreu et al.1 main ideas will be used in the definition of our new 
Chen and Wul d2001bh.lChan et alJ ^OO^.lNikolov^ ter, we introduce our filter. In section H we provide 

visual examples and numerical results that demonstrate 
our method's soundness. Sectionals a brief conclusion. 



where B is the set of pixels contaminated by impulse 
noise, P(B) = p is the impulse probability (the propor- 
tion of the occurrence of the impulse noise), n(x) are 
independent random variables uniformly distributed on 
some interval [c, d] . The impulse noise is characterized 
by replacing a portion of an image's pixel values with 
random values, leaving the remaining one unchanged. 
Such a noise can be introduced due to transmission er- 
rors, malfunctioning pixel elements in the camera sen- 
sors, faulty memory locations, and timing errors in analog- 
to-digital conversion. Recently, some important meth- 
ods have bee n proposed to remove t he impulse noise, see 



ocls nave been proposed to remove t ne imp i 
for ex a mple: Hwang and Haddadl ( 19951). 



Wenbini (|200.il JDong et al.l(|2007h . and 



Yu et all 



To remove a mixture of the Gaussian and impulse 
noises which defined by 



Y(x) = 



n(x), if x G B, 

f (x) + e(x), if i G I\B, 



(3) 



2 Algorithms 

2.1 Optimal Weights Filter 

G I and a given h > 0, the square 



For any pixel Xq 
window of pixels 



U 



x ,h 



{xel 



xo\ 



< h} 



(4) 



the above mentioned methods are not effective: the 
Gaussian noise removal methods cannot adequately re- 
move impulse noise, for they interpret the impulse noise 
pixel as edges to be preserved; when impulse removal 
methods are applied to an image corrupted with the 
Gaussian noise, such filters, i n practice, leave gra iny, 
visually disappointing results. iGarnett et al.l (|2005l) in- 
troduced a new local image statistic called Rank Or- 
dered Absolute Difference (ROAD) to identify the im- 
pulse noisy pixels and incorporated it into a filter de- 
signed to remove the additive Gaussian noise. As a re- 
sult they have obtained a trilateral filter capable to re- 
move mixed Gaussian and impulse noise. This method 
also performs well for removing the single impulse noise. 
For other developpem e nts in this dire c tion we refer to 

Boubouhs et all & . lLi1^D M . IXiao et al.l il2m7R ,1( ' r tobe distinguished from the search window U xo , h . 

The size of the patch window ~V x>r , is the positive in- 
teger m = (2n + l) 2 = card V XOi ^. The vector Y^ j?? = 
(Y (y))j /gV formed by the values of the observed noisy 



will be called search window at xq, where h is a positive 
integer. The size of the square search window \J XOl h is 
the positive integer number M = (2h+l) 2 = card \J Xo ,h- 
For any pixel x G U xo ,h and a given integer rj > a 
second square window of pixels 

Vj,, = {y G I : \\y - alloc < v} (5) 
will be called for short a patch window at x in or- 



Luisier et all |201ll) . 



In this paper, we propose a new filter that we call 
Optimal Weights Mixed Filter (OWMF). The idea of 
our method come from the combination of the ROAD 

and the Optimal Weights 



statistic of Garnett et al 



Filter in I Jin et al.1 (|2011l) . In the paper, we modify the 
Rank-Ordered Absolute Differences (ROAD) to Rank- 
Ordered Absolute Differences of mixture of Gaussian 
and impulse noises (ROADG). It will be more effective 
to detect impulse noise when the impulse is mixed with 



image Y at pixels in the patch ~V XyTj will be called sim- 
ply data patch at x G \5 Xo ,h- For any xq G I and 
any x G V Xa ,h, a distance between the data patches 



- X,7] 



= (Y(y)) 



defined by 



and Y 



Xo,V 



= (Y(y)) 



d 2 (Y 



x, rii Y^o.r;) 



Y x , i] 



is 



(6) 
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where 



£ {Y{T x y)-Y{y)f 



and T x is the translation mapping: T x y = x + (y — xq). 
As 

Y(T x y) - Y(y) = f(T x y) - f(y) + e(T x y) - e(y) 
we have 

E(Y(T x y) - Y{y)f = (f(T x y) ~ f{y)f + 2a 2 . 

If we use the approximation 

(f(T x y) - f{y)f « (f(x) - f(x )) 2 = P ), Xo {x) 

and the law of large numbers, it seems reasonable that 

pl Xo (x)*d 2 (Y x , v -Y Xo , 71 )-2a 2 . 



As shown in ([Jin et all 120111 ). in practice, a much bet- 
ter denoising results are obtained by using the following 
approximation 



Pf,x (x) ~ Px {x) = (d(Y XyTI ~ Y Xo>v ) - V2cr) 



(7) 



The fact that p XQ (x) is a good estim ator of Pf Xn was 
justified by the convergen ce results in (IJin et all 120111) 
(cf. Theorems 3 and 4 of (jjin et all l201lK The Opti- 
mal Weights Filter is defined by 



E KtrC^m*) 



OWF(/)(X ) = 



i£U, 



0.'' 



E 



(8) 



(9) 



where K tr is the usual triangular kernel: 

KtI (t) = (l - \t\)+ , ten 1 . 

Remark 1 The bandwidth a > is the solution of 

Px {x){a ~ p X0 (x)) + = a 2 , 



and can be calculated as follows. We sort the set {p Xo (x) : 
x € U Xo .h} in the ascending order = p Xo (x\) < 
Px {x2) <■■ < Px (x M ) < Px ( x m+i) = +00, where 
M = cardtXro h- Let 



a k 



and 



0-2 + E %o( x i 
i=l 

fe 

E Px ( x i) 



1< fc < M, 



(10) 



fc* = max{l < k < M : a k > p Xo (x k )} 
= min{l < k < M : a k < p xo ( x k)} - 1, 



(11) 



with the convention that = 00 if /^(xfc) = and 
that min0 = M + 1. The solution can be expressed 
as a — a k * ; moreover, k* is the unique integer k G 
{1, ••■ , M} such that aft > p XQ {x k ) anda fe+ i < /^(xfe+i) 
if fe < Af . 



The proof of Remark \T\ can be found in (jjin et al 



2011). 



2.2 ROAD statistic and Trilateral Filter 



In ( Garnett et al. . 20051) . Garnett et al introduced the 
Rank-Ordered Absolute Differences (ROAD) statistic 
to detect points contaminated by impulse noise. For any 
pixel xq £ I and a given d > 0, we define the square 
window of pixels 

^ ,d={x-0<N\\x-x \\ oo <d}, 

where d is a positive integer. The square window will be 
called deleted neighborhood at xq . The ROAD statistic 
is defined by 



A 



x G 1, 



(12) 



where ri(xo) is the i-th. smallest term in the set {|y(x)- 
Y(x n )\ 



x G f%LA and 2 < K < card 



x ,d' 



In 



( Garnett et all 120051 ) it is advised to use d = 1 and 
K = 4. Note that if Xq is an impulse noisy point, the 
value of ROAD(xq) is large; otherwise it is small. 

Following [Garnett et all (|2005l ) and lLi et all (|2010l ) 
the authors define the "joint impulsively" Jj (xq, x) be- 
tween xq and x as: 



J 1 (xq , x) = exp 



(ROAD(x Q ) + fiOAJ(x)) ; 



(13) 



where the function Jj(xq,x) assumes values in [0,1] 
and the parameter <jj controls the shape of the func- 
tion Jj{xq,x). If xq or x is an impulse noisy point, 
then the value of ROAD{xq) or ROAD(x) is large and 
Ji(xo,x) ~ 0; otherwise, the value of ROAD(xq) and 
ROAD(x) are small and Jj (xp,x) ~ 1. The definition 
of the trilateral filter (cf. ( Garnett et all [2005] )) is given 
by 



TnF(v)(x ) = 



4 
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where 

w(x) = ws(x)w R (x) Jl( - X0 ^w I (x) 1 ~ Jl( - x ''- x \ 

i*-»nl 2 

ws(x) = e 2 "s , 
^ii(^) = e 



The smoothing kernels k used in the simulations are 
the Gaussian kernel 



K g (y) = exp 



(y(x)-v-(x ))^ 

-ROAD(x) 2 
^ 



(17) 



wi(x) ■■ 

This filter has been shown to be very efficient in remov- 
ing a mixed noise composed of a Gaussian and random 
impulse noise. 



2.3 Optimal Weights Mixed Filter 



where h g is the bandwidth parameter and the following 
kernel: for y £ U XOtV , 



ko (y) 



E 



i 



<(ij) 



(2k + iy 



(18) 



The ROAD statistic (cf. lGarnett et al.l (|2005h ) provides 
a effective measure to detection the pixel contaminated 
by impulse. In this paper, we take into account the char- 
acter of Gaussian noise, and modify the method ROAD 
to adapt to the mixture of impulse and Gaussian noises. 
Then the equation (IT2l becomes 



, x £ I, (14) 



ROADG(x ) = (^J2 ri ^~ a 



where a is the standard deviation of the added Gaussian 
noise, ri(xo) is the i-th. smallest term in the set {|Y(a:) — 
Y(x )\ : x € ^° 0jrf }, and 2 < K < card fi° od . Let 



J(x, H) = exp 



ROADG(x) 2 
IP 



(15) 



be a weight to estimate whether the point is impulse 
one, where the parameter H controls the shape of the 
function. If the pixel x is an impulse point then ROADG(x) 
is large and J(x, H) w 0; otherwise ROADG(x) « 
and J(x, H) sa 1. 



if \\y — xq\\oo = j for some j € {0, 1, • • • , 77}. In the 
simulations presented below we use the kernel n — kq. 

Now, we define a new filter, called Optimal Weights 
Mixed Filter (OWMF), by 

M x °> = d; — ; ; — 7s~, — ixu — ' ( 19 ) 

2/£U XQ>fe 

where the bandwidth aj > can be calculated as in 
Remark[T](with p Xo (x) and a replaced by pj tK ,x ( x ) and 
aj respectively) and H2 is a parameter. Notice that Hi 
and H2 (which is used in the definition of p,j^ K . Xo (x)) 
may take different values. 

To explain the new algorithm (fT9)) . note that the 
function J(x,H2) acts as a filter of the points contam- 
inated by the impulse noise. In fact, if x is an impulse 
noisy point, then J(x, H2) ~ 0. When the impulse noisy 
points arc filtered, the remaining part of the image is 
treated as a image distorted by solely the Gaussian 
noise. So, in the new filter , the basic idea is to apply 
the OWF (jjin et all 120111 ) by giving nearly weights 
to impulse noisy points. 



Now, we modify the Optimal Weights Filter (jjin et al 



2011) in order to treat the mixture of impulse and Gaus- 
sian noises. Similar to ©, we define the impulse detec- 
tion distance by 



where 



I 0^x,ri ~^x ,ri)\\ 



J,K 



Ejl' 



K(y') 



= J2 <T x y)J(T x y,Hi)J(y,Hi)(Y(T x y)-Y(y)) 2 , 

and k are some weights defined on V Ko r; . The corre- 
sponding estimate of brightness variation Pf, Xo (x) is 
given by 



Pj.k,x {x) = (d JiK (Y x . v ,Y X(uV ) - V2c^j 



(16) 



3 Simulation and comparisons 

The performance of a filter / is measured by the usual 
Peak Signal-to- Noise Ratio (PSNR) in decibels (db) de- 
fined by 

P5^=101og 10 — , 

MSE = ^diT.(f^-f^ 2 > 

x£l 

where / is the original image. 

In the simulations, to avoid the undesirable border 
effects in our simulations, we mirror the image outside 
the image limits. In more detail, we extend the image 
outside the image limits symmetrically with respect to 
the border. At the corners, the image is extended sym- 
metrically with respect to the corner pixels. 

Algorithm : Optimal Weights Mixed Filter 
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For each x € I 



compute ROADG(x) = ( -g> T^f=i r i( x ) ~ a 

, j, tt \ ? ROA£>G(x) 

compute J (a;, Hi) = exp I v ; 

compute J(x,H2) — exp I H $ v ; 

Repeat for each x$ e I 

give an initial value of a: a = 1 (it can be an 
arbitrary positive number). 

compute {pj tK , Xo (x) : x <E U l0> 4 by (JTSJ) 
/ compute the bandwidth a at Xo 

reorder {pj, K , XQ (x) : x € U Xo ,h} as increasing 
sequence, say 

Pj,k,x (xi) < pj, K , Xo (x 2 ) < ■ ■ ■ < Pj, k , xo (xm) 

loop from k = 1 to M 

(xi) > 

lf E?-iP J .-. (^) " pJ ^ Xk) then 
- = ^ 2 +lXiP./.^o(*.) 

else quit loop 
else continue loop 
end loop 

/compute the estimated weights w at xq 



20% 



:30% 



40% 



50% 



compute w(xi) 

/compute the filter fh at Xq 

compute fh(x ) = J2 Xl eUx . h w(xi)Y(xi). 

In our simulations, the parameters can be choose as 
follows: 

d = 2, 
K = 12, 
M = 13 x 13, 
m — 15 x 15, 

*- + 

H 2 = 27- 20p. 



+ (a-10)+(0.5-p), 



In (jGarnett et al.1 . 120051 ) it is suggested to take d = 1 



and K = 4. For low and moderate levels of noise (p < 
25%), one iteration is sufficient and usually provides the 
best results; for high levels of noise (p > 25%), apply- 
ing two to five iterations provides better results. In our 
simulations, we found that a few spots of unremoved 
impulses often remain if we choose d — 1 and K = 4. 
This happens because impulses sometimes "clump" to- 
gether, and the 3x3 detection window is too small to 
identify all the impulse noise points. Consequently, we 
select parameters d = 2 and K = 12 of detection win- 
dows for all levels of impulse noise. Figure [T] shows the 
comparison results between the restored images, with 




PSN R=36 .03dbPSNR=33.65dbPSNR=30.22dbPSNR=27.35db 




Hill 

PSNR=35.50dbPSNR=33.92dbPSNR=32.19dbPSNR=30.09db 



Fig. 1 The first row gives the levels of added impulse noise. 
The second and the third rows display images restored with 
detection window 3x3, and their PSNR values. The forth and 
fifth rows show images restored with detection window 5x5 
and their PSNR values . The 1, 2, 3 and 4 columns give the 
restored images which have been contaminated by an impulse 
noise with p = 20%, 30%, 40%, and 50% respectively. 




Fig. 2 The first and second 100 x 100 images are parts of 
images restored with detection windows 3x3 and 5x5 re- 
spectively; the original image has been contaminated by an 
impulse noise with p = 20%. 



detection window 3x3 and with detection window 5x5, 
which have been added an impulse noise with p = 20%, 
30%, 40%, and 50% respectively. When p = 30%, 40% 
and 50%, we can see clearly some impulse spots in the 
restored images with detection window 3x3, while the 
visual quality of the restored images with detection win- 
dow 5 x 5 is very good, without impulse spots. In the 
case where p = 20%, impulse spots of the restored im- 
age with detection window 3x3 are not obvious, and 
the PSNR value is a little better than that with de- 
tection window 5x5, whereas Figure [5] shows that the 
first image has two clumpy impulse spots and the visual 
quality is not good enough. Consequently, we prefer de- 
tection window 5x5 for all levels impulse noise. 

We use the kernel ko for computing the estimated 
brightness variation function pj, K , X01 which corresponds 
to the Optimal Weights Mixed Filter as defined in Sec- 
tion [231 The parameters m and M have been fixed to 
m = 25 x 25 and M = 13 x 13. In Figure^ the images 
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Table 1 Performance of denoising algorithms when applied to test Gaussian white noisy images. 



Images 
Sizes 



Lena 
512 X 512 



Barbara 
512 X 512 



Boat 
512 X 512 



House 
256 X 256 



Method ~~ 
Our method 
M = 13 x 13 
25 x 25 



PSNR 



PSNR 



PSNR 



PSNR 



Buadc s et al. ('20051 
15 



20 



25 



33.75db 



31.81db 



31.02db 



32.72db 

Katkovnik ct al. (2004) 
Foi et al. (2004) 
Roth and Black (2009) 
Hirakawa and Parks (2006) 
Kcrvrann and Boulaneer (2008) 
Jin et al. /Mill 
Hammond and Simoncclli (2008) 
Aharon et al. (2006) ~ 

13abov et al (20o7T 

Our method 
M = 13 X 13 
m = 27 x 27 



31.67db 
32.18db 
32.72db 
33.29db 
33.97db 
33.70db 
33.93db 
34.04db 
33.71db 
34.27db 



30.39db 
29.20db 
29.61db 
30.16db 
32.55db 
31.80db 
32.31db 
32.25db 
32.41db 
33.00db 



33.82db 
30.46db 
30.93db 
31.27db 
31.59db 
31.44db 
31.64db 
31.72db 
31.77db 
32.14db 



32.42db 



30.40db 



29.62db 



Buadcs ct al. (2005) 
Katkovnik~ al. (2004') 
Foi et al. (2004) 
Roth and Black (2009^ 
Hirakawa and Parks (2006) 
Kervrann and Boulangcr (2008) 
Jin et al. (2011) 

Hammond and Simoncclli (2008) 
Aharon et al. (2006) 
Uabov et aT*(2007) 



Our method 
M = 13 X 13 
27 x 27 



31.51db 
30.74db 
31.43db 
31.89db 
32.69db 
32.64db 
32.68db 
32.81db 
32.39db 
33.05db 



30.38db 
27.38db 
27.90db 
28.28db 
31.06db 
30.37db 
31.04db 
30.76db 
30.84db 
31.78db 



29.32db 
29.03db 
39.61db 
29.86db 
30.25db 
30.12db 
30.30db 
30.41db 
30.39db 
30.88db 



31.40db 



29.20db 



28.56db 



Buadcs ct al. (2005) 

Katk ovnik et aT (2004) 

Foi et al. I'iMi] 

Roth and Black (,2009) 

Hirakawa and Parks (2006) 

Kcrvrann and Boulangcr (2008) 

Jin et al. (2011) 

Hammond and Simoncclli (2008) 

Aharon et al. (2006) ~ 

Uabov et al . (2007) 



30.36db 
29.66db 
30.43db 
30.57db 
31.69db 
31.73db 
31.59db 
31.83db 
31.36db 
32.08db 



29.19db 
26.05db 
26.62db 
26.84db 
29.89db 
29.24db 
29.92db 
29.58db 
29.58db 
30.72db 



28.38db 
27.93db 
28.60db 
28.57db 
29.21db 
29.20db 
29.16db 
29.40db 
29.32db 
29.91db 



33.82db 



32.62db 
33.18db 
33.55db 
33.82db 
34.08db 
34.09db 
33.72db 
34.25db 
34.94db 



32.71db 



32.51db 
31.24db 
31.84db 
32.29db 
32.58db 
32.90db 
32.83db 
32.52db 
33.10db 
33.77db 



31.61db 



31.16db 
30.12db 
30.75db 
31.05db 
31.60db 
32.22db 
31.95db 
31.54db 
32.07db 
32.86db 



Table 2 Performance of denoising algorithms when applied to test impulse noisy images. 



Images 



Baboon 



Bridge 



Lena 



Pentagon 



2% 



"20%" 



"40%" 



"20% 



"20%" 



"40%" 



"20% 



40%" 



Method 



PSNR PSNR 



PSNR PSNR 



PSNR PSNR 



PSNR PSNR 



Our method 
M = 13 X 13 
25 x 25 



24.81db 22.12db 



27.84db 24.91db 



35.50db 32.19db 



30.91db 28.34db 



Sun and Ncuvo (1994) 
Abreu et aTTTbwFi 
Wan^inc^hanE; (1999) 
Ohen et al. ( 1999) 
C"hen and WuT^OOlbl 
Chen and ^u (2001a ) 
(Jrnoievic et alM^OW) 
Wenbin (26T5T 
Garnctt ct al. (2005) 
Chan ^t"al~72004) 
Dong et al" (2007) 
Vu et al. (2*008) 



23.67db 
23.81db 
23.43db 
23.73db 
24.02db 
24.17db 
23.78db 
24.18db 
24.18db 
23.97db 
24.49db 
24.86db 



20.85db 
21.49db 
21.07db 
21.38db 
21.52db 
21.58db 
21.56db 
21.41db 
21.60db 
21.62db 
21.92db 
22.06db 



26.26db 
26.56db 
26.33db 
26.52db 
27.27db 
27.08db 
26.90db 
27.05db 
27.60db 
27.31db 
27.86db 
28.06db 



22.66db 
23.80db 
22.75db 
22.89db 
23.55db 
23.23db 
23.83db 
23.88db 
24.01db 
24.60db 
24.79db 
24.97db 



32.93db 
35.71db 
35.09db 
34.21db 
35.44db 
36.07db 
36.50db 
36.90db 
36.70db 
36.57db 
37.45db 
36.18db 



27.90db 
29.85db 
28.92db 
28.30db 
29.26db 
28.79db 
31.41db 
30.25db 
31.12db 
32.21db 
32.76db 
32.03db 



29.34db 
30.38db 
29.18db 
29.29db 
30.34db 
30.23db 
30.11db 
30.42db 
30.33db 
30.03db 
30.73db 



26.26db 
27.27db 
26.19db 
26.29db 
27.04db 
26.84db 
27.33db 
26.93db 
27.14db 
27.35db 
27.73db 



in the third row show that the noise is reduced in a 
natural manner and significant geometric features, fine 
textures, and original contrasts are visually well recov- 
ered with no undesirable artifacts. To better appreciate 
the accuracy of the restoration process, the images of 
square errors (the square of the difference between the 
original image and the recovered image) are shown in 



the fifth row of Figure |3l where the dark values corre- 
spond to a high-confidence estimate. As expected, pix- 
els with a low level of confidence are located in the 
neighborhood of image discontinuities. For comparison, 
we show the images denoised by the trilateral filter TriF 
(see the images in the second row of Figure [3] ) and 
their square errors (see the images in the forth row of 
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Table 3 Comparison between the TriF dGarnett et al.l . 120051 ). and OWMF (Our method) in removing mixed noise. 



Gaussian Moisc 


Image 


Method 


p = 0.2 p = 0.3 p = 0.4 p = 0.5 


sigma=10 


Lena 


Garnctt et al. (2005) 
Our method 


31.48db 29.87db 28.57db 27.31db 
33.18db 32.05db 30.90db 29.52db 


Bridge 


Garnctt et al. (2005) 
Our method 


25.82db 24.92db 23.79db 22.28db 
26.42db 25.19db 24.08db 23.08db 


Boat 


Garnctt et al. (2005) 
Our method 


28.61db 27.54db 26.22db 24.74db 
29.57db 28.22db 27.05db 25.92db 


Barbara 


Garnctt ct al. (2005) 
Our method 


24.82db 24.00db 23.08db 22.33db 
28.47db 26.46db 24.83db 23.62db 


sigma=20 


Lena 


Garnctt ct al. (2005) 
Our method 


28.85db 28.02db 27.10db 25.68db 
30.87db 30.09db 29.19db 28.14db 


Bridge 


Garnctt ct al. (2005) 
Our method 


23.56db 23.01db 22.47db 21.72db 
24.70db 23.97db 23.21db 22.45db 


Boat 


Garnctt ct al. (2005) 
Our method 


26.18db 25.46db 24.75db 23.79db 
27.79db 26.93db 25.97db 25.08db 


Barbara 


Garnctt ct al. (2005) 
Our method 


23.35db 22.95db 22.53db 21.84db 
27.50db 25.95db 24.43db 23.33db 


sigma=30 


Lena 


Garnctt ct al. (2005) 
Our method 


27.26db 26.57db 25.58db 23.99db 
29.12db 28.49db 27.76db 26.75db 


Bridge 


Garnctt ct al. (2005) 
Our method 


22.88db 22.42db 21.87db 20.98db 
23.56db 23.02db 22.49db 21.86db 


Boat 


Garnctt ct al. (2005) 
Our method 


25.11db 24.55db 23.80db 22.ti2db 
26.41db 25.79db 25.08db 24.26db 


Barbara 


Garnctt ct al. (2005) 
Our method 


22.82db 22.46db 21.94db 21.10db 
25.98db 24.81db 23.72db 22.81db 



Figure [3j . We can see clearly that the images of square 
errors of our methods are darker than that of Trif, so 
our method provides significant improvement over. The 
overall visual impression and the numerical results are 
improved using our algorithm. 

For comparison, we consider follows three cases: pure 
Gaussian noise, pure impulse noise and the mixture 
of Gaussian and impulse noises. In the case of pure 
Gaussian white noise, we have done simulation on a 
commonly- used set of images ("Lena" , "Barbara" , "Boat" 



(Il996h IWang and Zhanel fll999h IChen et all (ll999l).IChen and Wu 
d2001bl) IChen and Wul (|2001ah . lCrnoievic et al.l ( 2004h . 
IWenbinl(l2005ir etc. Our app roach work as well as lDong et al 




and "House") available at |http: / / decsai.ugr.es/~javier / 



denoise/test_images/ and the comparison with several 
filters is given in Table [TJ The PSNR values show that 
our app roach work as well as those s ophisticated meth- 



(|2007T) and lYu et al.l (|2008[ ) , when our approach pro- 
duces the best PSNR values in the cases of " Baboon" 
40%) and "Pentagon" (20% and 40%), while I Yu et al 
has the best results in the case o f "Baboon" 
and "Bridge" (20% and 40%), and bong et"aH 
(ROLD-EPR) wins in the case of "Lena" (20% 
and 40%). Finally , we s how the comparison between 
the iGarnett et al.l (|2005l) and our method with the set 
of images ("Lena", "Bridge", "Boat" and "Barbara"). 
From Table [HJ it is clear that our method pr ovides sig- 



ods, li ke Hirakawa and Parks! (l2006l).lKervrann and Boulangeffi ^ 



(2008), H ammond and Simonce lli (2008) and lAharon et al 



mficant improvement over the algorithm Trif (|Garnett et al 



( 2006) , and is better than the filters proposed in 



(2005), Salmon and Lc 



ennee 



Buades et a' 



(2009), Katkovnik et al 



Roth and Blackl <|2009h 



Conclusion 



20041) . iFoi et al.l (|2004l) and 

The proposed approach gives a quality of denoising 
which is co mpetitive with t hat of the best recent method 



BM3D (see babov et al.l |2007t )K These methods can 



only deal wits pure Gaussian noise, while our method 
can not only cope with the Gaussian noise, but also 
remove the impulse noise and the mixture of Gaussian 
and pure impulse noises. For the pure impulse noise, our 
method is also competitive to the sophisticated method. 
In order to compare the others methods, we choose a 
commonly set of images ("Baboon" , "Badge" , "Lena " 



A new image denoising filter, which deal with the mix- 
ture of Gaussian and impulse noises model based on 
weights optimization and the modified Rank-Ordered 
Absolute Differences statistic, is proposed. The imple- 
mentation of the filter is straightforward. Our work 
leads to the following conclusions. 

1. The modified Rank-Ordered Absolute Differences 
statistic, used in the new filter, detects effectively 
the impulse noise in the case of mixture of Gaussian 
and impulse noises. This statistic is well ada pted for 



use w ith the Weights Optimization Filter of (|Jin et al 



2011) 



and "Pentagon") which is taken in (|Dong et all 12007') . 
Table [5] lists the restoration results of well- know differ- 
ent algorithm. It is clear that our method provides sig- 2. The proposed filter is proven by simulations to be 
nificant improvement over Sun and Neuvol ( 1994 ) , Abreu et al. very efficient for removing both a mixture of impulse 
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Qiyu Jin et al. 




Fig. 3 The first row is the images with different levels of 
mixture of impulse and Gaussian noise. The second and the 
third row display the images restored by TriF and OWMF 
respectively. The forth and the fifth row give square error 
image with TriF and OWMF respectively. 

and Gaussian noises, and the pure impulse or pure 
Gaussian noise. 
3. Our numerical results demonstrate that the new fil- 
ter outperforms the known filters. 
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